Orbital Dimerization and Dynamic Jahn- Teller Effect in NaTiSi2 06 
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, We study with Raman scattering teclinique two types of phase transitions in the pyroxene corn- 

ed ' pounds NaMSi206 (with M=Ti, V, and Cr). In the quasi one-dimensional S=l/2 system NaTiSi206 

we observe anomalous high-temperature phonon broadening and large changes of the phonon ener- 
gies and line-widths across the phase transition at 210 K. The phonon anomalies originate from an 
■ orbital order-disorder phase transition and these results -combined with theoretical considerations- 

indicate that the high temperature dynamical Jahn- Teller phase of NaTiSi2 06 exhibits a sponta- 
neous breaking of translational symmetry into a dimerized, Jahn- Teller distorted, orbital ordered 
^ . state under the formation of spin valence bonds. In S=l NaVSi206 orbital degrees of freedom are 

strongly suppressed and the magnetic excitations are well described within a Heisenberg model, 
indicating that at Tjv=19K this system orders antiferromagnetically. 
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I , Introduction. Electrons in strongly correlated transition-metal compounds can be regarded as having separate 
' ^ ■ spin, charge and orbital degrees of freedom. It is the interplay between these, combined with their coupling to the 
^ lattice, that gives rise to a wealth of possible spin, charge and orbital orderings, as observed for instance in many 
Q Colossal Magneto- Resistance manganites 1], cuprates titanates (e.g. LaTiOa Q) and vanadates (e.g. V2O3 0| 
I 1[ and LiVOa 15|). 

Systems with orbital degeneracy are particularly interesting because orbitals couple to the lattice via the cooperative 
' Jahn- Teller (JT) effect on one hand, and via superexchange interactions to the electronic spin on the other hand 0]- 
I , Therefore at an orbital ordering phase transition the magnetic susceptibility and phonon properties will be affected 
at the same time. Experimentally, however, such an interrelation is seldom found: in general the dominant JT 
T— I orbital-lattice coupling obscures the more subtle effects due to the superexchange. 

Pyroxenes are large family of compounds (AMB2O6; A=alkali-metal, M=transition- metal, and B=Si, Ge) which 
^""^ ] structure consists of isolated quasi one-dimensional chains of edge-sharing MOe octahedra, linked together by the 
. corner-sharing BO4 tetrahedra, Fig. 1. Particularly, the compounds with active t2g orbitals (M=Ti,V, and Cr) 
are expected to have electronic interactions governed by both orbital degeneracy and anisotropy. Interestingly, the 
sodium-silicon system with titanium, NaTiSi206, is rather different from other pyroxenes as it lacks low-temperature 
AF order, and shows signs of the opening of a spin gap instead 0,0 Ti'^+ corresponds to S=l/2, and since 
all Ti ions are in equivalent crystallographic positionSjthis material is a prime candidate to show a spin-Peierls (SP) 
phase transition. Indeed, its magnetic susceptibility y\ sharply decreases below 210 K, indicating a transition to a 
^ _ spin-singlet state, which could in principle be of the SP type. But already the high-temperature magnetic susceptibility 
O ■ data does not support such a SP scenario, since the phase transition occurs at a temperature that is higher then the 
maximum point of the Bonner-Fisher curve, which implies that the transition cannot solely be driven by magnetic 
fluctuations 6]. Based on structural analysis, Isobe et.al. anticipated that orbital dimerization might be responsible 
for such behavior 

^ , In this Letter, we report a Raman scattering study on a family of pyroxene compounds from which we conclude 
■ that orbital degrees of freedom play decisive role in physical properties of NaTiSi206, and that NaTiSi206 undergoes 
an orbital order phase transition at Too=210 K. At high temperatures the fluctuations of the orbital degrees of 
freedom lead to a dynamic Jahn- Teller phase with anomalous phonon broadening and remnant antiferromagnetic 
(AF) spin fluctuations. We find that the dramatic drop of the magnetic susceptibility below Too is accompanied 
with a structural change, just as one would expect for a canonical orbital ordering transition. These results, combined 
with the microscopic orbital-spin model that we derive, suggest that the quasi one-dimensional dynamical Jahn- Teller 
phase of NaTiSi206 exhibits a spontaneous breaking of the translational symmetry into an dimerized orbital ordered 
state with a spin gap due to the formation of spin valence bonds. The paper is organized as follows: first we present 
the Raman spectra and underline important observations, then we discuss the electronic structure of pyroxenes and 
derive a microscopic model. Finally, we demonstrate that all observed effects (including susceptibility data) can be 



2 




accounted by our model. 

ExperimentMigh-quality powder samples of NaMSi206 (M=Ti, V, and Cr) were prepared by a solid-state reaction 
of mixtures with an appropriate molar ratio of Na2MSi40ii, M and M02- Details of the sample preparation are 
published elsewhere |^]. Raman spectra are measured in the backscattering configuration, using 514.5 nm line of an 
Ar-ion laser, micro- Raman system with a DILOR triple monochromator, and a liquid nitrogen cooled charge-coupled 
device detector. For the low-temperature measurements we used the Oxford continuous-flow cryostat with a 0.5 mm 
thick window. The focusing of the laser beam is realized with a long distance (10 mm focal length) microscope 
objective (magnification 50 x). 

The pyroxenes crystallize in a monoclinic unit cell with the space group C2/c 0. The unit cell consists of four 
formula units (Z=4) with 40 atoms in all. The site symmetry of Na, M, Si, Oi, O2 and O3 atoms are (4e), (4e), (8f), 
(8f), (8f) and (8f), respectively. Thus, the factor-group-analysis (FGA) yields: (Na, M)(C2) T ^ Ag + A,, + 2Big + 
2Bu; (Si, Oi, O2, 03)(Ci) F = 3Ag + 3A„ -I- 3Bg + SB^. Summarizing these representations and subtracting the 
acoustic modes {Au + 2Bu), we obtain the following irreducible representations of NaMSi2 06 vibrational modes: 

^NaMSi20e = 14Ag(x2;, yy, zz, xz) + \'oBg{xy, yz) 
+ m„(E||y)-f 14B„(E||x,E||z) 

The unpolarized Raman spectra of NaMSi206 (M=Ti, V, Cr) are shown in Fig. 2. At the room temperature we find 
around 30 phonon modes as predicted by FGA. The Raman spectra of different compounds are similar as expected 
for isostructural materials, and the phonon modes can be crudely grouped into two energy regions. The modes 
in the spectral range below 500 cm^-'^ originate from the bond bending vibrations, whereas the higher frequency 
modes originate from the stretching vibrations. The highest energy modes are mainly due to the non-bridging Si-0 
ion vibrations, because of the shortest Si-0 tetrahedral bonds. The modes of NaMSi206 at 1055/1032 cm"'^ (Cr), 
1042/1025 cm~^ (V), and 1042/1025 cm^^ (Ti), represent Si-02 (see Fig. 2) antisymmetric/symmetric bond stretching 
vibrations, respectively. Their frequency difference scales as R~'^ in a full accordance with difference between Si-02 
bond lengths in these materials. Similarly, the modes at 990/967 cm~^ (Cr), 972/954 cm~^ (V), and 965 cm~^(Ti), we 
assign as antisymmetric/symmetric pairs of Si-Oi bond stretching modes. This simplified mode assignment becomes 
inapplicable at lower frequencies due to more complicated normal coordinates of corresponding vibrations. However, 
besides similarities, we also observe two very important effects. First, most of the phonon line widths are dramatically 
increasing in the components with smaller spin values {V^^ — > 5 = l,Cr'^+ S = 3/2). Second, the spectra of 
NaTiSi2 06, due to the large phonon broadening, show effectively less phonon modes then expected by FGA, and 
observed in other pyroxenes. The latter effect may also be regarded as a consequence of a high temperature "higher 



3 



300 K NaMSi,0^ 




200 400 600 ^ 800 1000 

Raman shift (cm" ) 



FIG. 2: Raman spectra of NaMSi206. 

symmetry lattice state" of NaTiSi2 0e. 

By lowering the temperature we find a dramatic change in the phonon Raman spectra of NaTiSi206 around 210 K, 
see Fig. 3a, b. The modes at 221 and 209 cm~^ exhibit the anti-crossing behavior, see inset Fig. 3a, which indicates 
that these two modes belong to the same symmetry class below Tc- A similar elTect is observed for all other such 
pair of modes (one more pair is also visible in Fig. 3a at about 180 and 190 cm~^) in the spectra, which suggests 
that in the low-T phase of NaTiSi206 all Raman active modes are of the same (Ag) symmetry type. Furthermore, 
our reflectance measurements (will be published separately, Ref. 0) show that the infrared active phonon modes 
can be distinguished from the Raman modes, since the center of inversion remains to be the symmetry element in 
the low-T phase. Accordingly, we uniquely determine the space group of the low-T phase of NaTiSi206 to be PI. 
This space group requires the translation symmetry breaking of the TiOe chain (the tentative distortion pattern of the 
low-T phase is shown in Fig. 1) which is in agreement with preliminary X-ray diffraction and neutron scattering data 

m 

Moreover, we show in details the temperature dependence of the structure around 970 cm ^, Fig. 3b. The mode 
at about 946 cm~^ softens by about 10 cm~^, while mode at 966 cm~^ "splits", and hardens by 25 cm~^, see inset 
in Fig. 3b. The full width at half maximum (FWHM) of the 946 cm~^ (circles) and 966 cm~^ (squares) phonons 
(the FWHM is estimated from Lorenzian-fit) increase up to the maximum value at about 210 K, and then decreases 
to the saturation value which is much smaller then the T=300 K value, see right inset of Fig. 3b. In fact, all 
Raman-active phonons exhibit anomalies at the temperature which coincides with Tc obtained from susceptibility 
measurements |6| indicating that the magnetic ordering is accompanied with the structural phase transition just as 
one would expect for an canonical orbital-ordering transition. Besides, behavior of FWHM implies that the bond 
fluctuations are considerably larger in the high-T then in low-T p/iase(proximity of the structural phase transition 
induces the largest fluctuations, producing the maximum FWHM at Too)- Due to the JT effect bond fluctuations 
emphasize the strong orbital character of disorder above the phase transition temperature 12, 13] . 

Contrary to S=l/2 compound, in the low temperature Raman spectra of NaVSi206, see Fig. 3c, we did not 
observe the phonon anomalies. Instead, at low temperatures, broad asymmetric features, typical for the two-magnon 
excitations in Heisenberg AF {H = X^i j Ji.j^i ' ^j^Jij > 0)i appear around 260 and 400 cm^^. According to the 
linear Anderson's approximation, differential cross section for the two-magnon (tm) Raman scattering is proportional 
to the two-magnon density of states (DOS), Itm ~ Sk'^('^ ~ 2w(k)). The dispersion relation is obtained assuming 
two-dimensional magnetic structure with different exchange interactions along the chains (J||), and perpendicular (Jj_) 
to the chains: w(k) — 25^/(^11 + J±)'^ — {J\\cos{ka) + J±cos{kb))'^ . The calculated two-magnon DOS is shown in top 
panel of Fig. 3c. Two singularities at energies of about 260 and 390 cm~^ are obtained for J|| = 85 cm~^ and Jj_ — 13 
cm~^ (120 K and 18 K, respectively), in a very good agreement with experiment. The maximum of the susceptibility 
curve m NaVSiaOe is at about 100 K [H (giving J,, - 80 K [l|) which is in good agreement with our estimation of 
the J||. Thus, we conclude that NaVSi20R exhibits an AF ordering. Similar AF phase has been already observed in 
other S=l chain compounds LiVGe206 and NaVGe206 0|- 



4 




0T(K)100 200 300 ^, 0T(K)100 200 300 0T(K)100 200 300 




700 750 800 850 900 950 1000 1050 1100 

Raman shift (cm"^) 




200 300 400 500 

Raman shift (cm"^) 



FIG. 3: Temperature dependence of the Raman spectra in a),b) NaTiSi206, and c) NaVSi206. Top panel in c) shows 
calculated two-magnon density of states. Insets: Temperature dependencies of the frequency and line width of various phonons 
in NaTiSi206. Vertical dashed line is an guide for the eye. 
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Pyroxenes: theoretical. Next we discuss the electronic structure. For a single octahedron the cubic crystal field 
splits the Ti,V or Cr 3d-states into low lying t2g states, with 1,2 or 3 electrons, respectively, and empty states of eg 
symmetry at higher energy. The low energy electronic properties are governed by the three-fold degenerate t2g states 
(the relevant states are \xy), \yz} and \zx)) [l^ . 

The Coulomb interaction U between electrons on the same transition metal atoms is large, so the exchange in- 
teractions can be determined by a second order perturbation expansion in the electron hopping parameters. The 
problem is further reduced by considering the symmetry allowed hopping paths in the chain geometry, schematically 
represented in Fig. 4. If we consider orbitals on two sites in the same xy plane, then only the hopping between \xy) 
orbitals is relevant (see Fig. 4a). For sites in the yz plane, the \yz) orbitals are relevant. Fig. 4b. In the present 
geometry (note that x, y, and z in Fig. 4 do not correspond to a, b, and c axis in Fig. 1) there are no transition metal 
atoms in the same chain that are also in the same xz plane. The \xz) orbitals are therefore non-bonding and can be 
considered inert on this level of approximation (a tightbinding parametrization |l7l| shows that other overlap integrals 
either vanish by symmetry or are more than factor five smaller), see Fig. 4c. For the S=l/2 Titanium system we 
then obtain the Hamiltonian [l8|. 

H^' =1 I J2 • [\ + T^TJ + t^{T^-+ T^)\ , (1) 

where we use the orbital operators T (Tf = 1/2 corresponds to an occupied \xy) orbital and T!f = —1/2 to an occupied 
\yz) orbital on site i), and i,j are neighboring sites. For the S=l Vanadium and S=3/2 Chromium system the orbital 
degree of freedom vanishes (T = 0) as the both of the active orbitals are occupied, so for the spin system we obtain 
a simple Heisenberg Hamiltonian iJ^/^'' =| J^/'-^'' | j ' ^i- This immediately explains the uniqueness of the Ti 
(S=l/2) system, and the antiferromagnetism in NaVSi206 (in fact in all S ^ 1/2 pyroxenes). 

Interpretation of experiments. The ground state of Hamiltonian JQ) , is clearly a ferro-orbital state, with spin-singlets 
on each bond, where the energy per dimer is — 3J'^'/4. The state with all \xy) occupied is degenerate with the state 
with all \yz) occupied, see Fig. 4. Those states do differ, however, because the dimerization pattern along the chain is 
shifted by one lattice spacing. At zero temperature the system is condensed in either one of these two dimerized orbital 
ordered states and the translation symmetry is broken. This explains the structural change at Tc ^Hj , the symmetry 
change and the energy shifts of the phonon excitations in the Raman spectra of NaTiSi2 06, and the observation of a 
large spin gap in the susceptibility measurements 

At high temperatures, due to strong orbital-lattice coupling and JT effect, the orbital fluctuations produce dra- 
matic effects. Above Too the JT distortion disappears: the high-T orbital disordered phase of NaTiSi206 may be 
regarded as an orbital fluctuating phase, and on larger time scales the TiOe octahedra appear to be undistorted. This 
causes the crystal to be effectively more symmetric, which is in agreement with the observation of less-then-expected 
phonon modes in the room temperature Raman spectra. As a consequence, orbital fluctuations induce large phonon 
broadenings (of course, the modes with Ti-0 bonds in their normal coordinates will be mostly affected), as they are 
indeed found in the room temperature Raman spectra of NaTiSi2 06, see Fig. 2. In that respect the high-T phase of 
NaTiSi2 06 resembles a dynamical J T phase (where the phonon broadening is a signature of the melted static lattice 
distortions (12.. .ISi]). The spin fluctuations above the orbital ordering temperature |i6|| can be easily understood in a 
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straightforward mean- field approximation of Hamiltonian where we can decouple orbitals and spins. If the long 
range order is absent in the orbital sector, we have for the expectation values {T^) = 0, and \{TIT^)\ < 1/4, so that 
the effective exchange constants -and therefore the fluctuations in the spin-sector- are still antiferromagnetic. Thus, 
the nature of the phase transition in NaTiSi206 can be established: it corresponds to an orbital order-disorder phase 
transition with appropriate concomitant magnetic and lattice changes. 

In the case of two electrons per site (NaVSi206), the orbital fluctuations are strongly suppressed due to inert 
property of \xz) orbitals, and the corresponding phonon broadenings are much smaller, see Fig. 2. NaCrSi206 has a 
fully polarized t2g core, no orbital degrees of freedom, and no anomalous phonon broadenings in the Raman spectra. 

Conclusions. We report a study of the two types of the phase transitions, observed in pyroxene family, by analyzing 
the Raman-active phonon and magnon excitations, and their temperature dependence. We find that spin S=l/2 
compound, NaTiSi206, exhibits a transition at about 210 K that we assign to be an orbital order-disorder phase 
transition. It originates from the instability of the high temperature orbital fluctuating - dynamical Jahn- Teller 
phase, towards a dimerized orbital ordered state, which is accompanied by a lattice distortion and by spin valence 
bond formation. The spin S=l system, NaVSi2 06, on the contrary, does not show this type of instability and orders as 
a Neel-type antiferromagnet below Tat = 19 K in agreement with our microscopic spin-orbital model for the pyroxenes. 
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